Wave-driven dynamo action in spherical MHD systems 
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Hydrodynamic and magnetohydrodynamic numerical studies of a mechanically forced two-vortex 
flow inside a sphere are reported. The simulations are performed in the intermediate regime between 
the laminar flow and developed turbulence where a hydrodynamic instability is found to generate 
internal waves with a characteristic m = 2 zonal wave number. It is shown that this time-periodic 
flow acts as a dynamo although snapshots of the flow as well as the mean flow are not dynamos. 
The magnetic fields' growth rate exhibits resonance effects depending on the wave frequency. Fur- 
thermore, a cyclic self-killing and self-recovering dynamo based on the relative alignment of the 
velocity and magnetic fields is presented. The phenomena are explained in terms of a mixing of 
non- orthogonal eigenstates of the time dependent linear operator of the magnetic induction equation. 
The potential relevance of this mechanism to dynamo experiments is discussed. 

PACS numbers: 91.25.Cw, 47.65.-d, 47.20.-k 



I. INTRODUCTION 



It is well-established theoretically that magnetic fields 
which emanate from the Earth and from celestial bodies 
are generated by flows of electrically conducting fluids 
in their interiors, a mechanism which is commonly re- 
ferred to as the magnetohydrodynamic (MHD) dynamo 
effect [H, In recent years, significant effort was de- 
voted to the experimental verification of dynamo theory 
in simply connected spherical [1, 0, Q and cylindrical 
[1, 01 impeller-driven flows of liquid sodium. While the 
latter setup has ultimately demonstrated a self-excited 
dynamo, it has turned out that turbulence strongly in- 
hibits the dynamo process in both experiments. A review 
on dynamo experiments is given in Ref. Detrimen- 
tal effects of turbulence on the excitation threshold of 
large-scale dynamos were confirmed by numerical simu- 
lations of MHD systems with well-defined mean-flows in 
periodic-box geometry lol. [Tol . [TTl [T2I as well as in spher- 
ically bounded geometry [l^, [iJ, [3 ■ 

In the present paper, we study numerically the regime 
between turbulence and the laminar state of an electri- 
cally conducting two- vortex flow inside a sphere which is 
similar to the flow realized in the Madison Dynamo Ex- 
periment 0, 13) Hi • The structure of the manuscript is as 
follows. The numerical model is presented in section [Til 
Section IIIII discusses purely hydrodynamic simulations 
which are performed in the sub-turbulent regime. It is 
found that a hydrodynamic instability leads to traveling 
internal waves, their dominant component being symmet- 
ric under rotations by tt. In section HVl we turn to MHD 
investigations of the flows introduced in the previous sec- 
tion. It is found that the presence of the waves which cor- 
respond to smoothly oscillating large-scale fluctuations 
supports the dynamo instability. Artiflcially changing 



the wave's frequency yields a resonance effect, i.e. the 
dynamo is found to operate most efficiently at certain 
frequencies whereas dynamo action ceases for frequencies 
which are too low or too high. Moreover, a nonlinear dy- 
namo is reported which undergoes a cycle of self-killing 
and self-recovering events. We find that this behavior is 
related to phase shifts between characteristic magnetic 
and hydrodynamic oscillations which are imposed during 
the transitions. These phase shifts translate to changes 
in the relative alignment of the velocity and magnetic 
fields. Section IVl concludes with a discussion and inter- 
pretation of the numerically obtained results. We inter- 
pret the MHD phenomena in the framework of a dynamo 
mechanism based on the mixing of non-orthogonal eigen- 
states of the time-dependent linear operator of the in- 
duction equation p^ . thus confirming its relevance to a 
self-consistent dynamo model. 



II. NUMERICAL MODEL 

The magnetic (B) and velocity (v) fields which de- 
scribe an incompressible electrically conducting fiuid are 
governed by the induction equation ([T]) and the Navier- 
Stokes equation ([2]), 
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along with the constraints V-B = V-v = 0. Here, A is the 
magnetic diffusivity, p is the pressure, j = /ig"^ VxB is the 
current density, v is the viscosity, and F is a forcing term. 
A constant mass density p = 1 is assumed. In nondimen- 
sional form, the problem is characterized by two control 
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parameters, the Reynolds number Re = LVu and the 
magnetic Reynolds number Rm = LVA"^, where L and 
V denote length and velocity scales which are characteris- 
tic of the system under consideration. We choose L = 1 
and V — Wrms, where Wrms = \/ (v^) is the spatial rms 
velocity. Here, the angle brackets denote averaging in 
space, and the overline denotes averaging in time which 
is performed during the quasi-stationary phase of the flow 
after all initial bifurcations. The characteristic timescale 
of the flow is given by the eddy turnover time = L/V, 
whereas the magnetic diffusion time t^- — Rmr^ is the 
timescale relevant to the magnetic field. 

We solve numerically the MHD equations ([I])-® in 
a sphere, using the parallel version of the Dynamo code 
[1J,[13- It employs the standard pseudo-spectral method 
based on a poloidal-toroidal decomposition of the vector 
fields in combination with spherical harmonic expansions. 
Outer boundary conditions arc the potential field solu- 
tion for the magnetic field, and, unless stated otherwise, 
the zero-slip condition for the velocity field. The forcing 
is designed to produce an s2t2 type of flow , which 
consists of two in the toroidal direction counter-rotating 
hemispherical cells with poloidal circulation in each cell, 
directed outwards to the poles and inwards in the equato- 
rial plane. It is given by an axisymmetric localized body 
force F kept constant in time, which reads in a cylindrical 
coordinate system (s, 0, z) 

Fs = 0, 

= e sgn(z)sV-3 sin ^ + 7, (3) 

TTS 

Fz — (1 — e) sgn(z) sin V 5. 

Assuming a sphere of radius r — 1, the driving is applied 
within 0.25 < \z\ < 0.55, s < = 0.29, with the param- 
eters kept constant at e = 0.1, 7 = 0.05, 6 = 0.3. For 
sufficiently large Rm, the resulting flow is a dynamo, the 
(£ = TO = 1) magnetic eigenmode growing fastest. The 
generation of this transverse dipole can be understood in 
the axisymmetric flow using a simple frozen-flux picture 
based on stretching, twisting, and folding of magnetic 
field lines 0. 

In all simulations, the spherical harmonic expansions 
are triangularly truncated at degree and order ^max = 
'^max = 20, while 160 uniformly spaced radial points 
are used. To verify the simulation results obtained at 
the aforementioned resolution we have repeated a nonlin- 
ear run (the computation of flow C which is introduced 
in the following section) using £max = Wmax = 30 and 
320 gridpoints in radial direction. The bifurcations to- 
wards flow C as well as the growthrate of the magnetic 
field in the MHD case were fully reproduced which in- 
dicates that converged solutions are obtained already at 
the coarser resolution. Moreover, power spectra of the 
velocity fields in terms of £ show a drop-off by three to 
four orders of magnitude between the dominant and the 
highest wavenumber. The magnetic field is resolved sim- 
ilarly well. Relevant spectra are shown and discussed in 
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FIG. 1: Time traces of Vims for four s2t2 flows driven by the 
body force F, at Re being moderately above the stable region 
Re < 65, illustrating bifurcation sequences. 

the following sections, cf. Figs. H] and [T^l After having 
introduced the numerical setup we now turn to an in- 
vestigation of the hydrodynamics of the flow before we 
consider full MHD cases. 



III. HYDRODYNAMIC STUDIES 

We flrst investigate the dynamics of the unmagnetized 
flow, i.e. Eq. ^ is integrated forward in time with 6 = 0. 
Starting from a fluid at rest, momentum is injected by 
the body force F, and the time integration is performed 
until a statistically stationary state is reached. This pro- 
cedure is repeated to scan over various decreasing vis- 
cosities, thus increasing Re. For Re < 65, the flow 
reaches a stationary axisymmetric (to — 0) state. Above 
this threshold, the flow becomes hydrodynamically un- 
stable. To exemplify subsequent temporal transitions. 
Fig. [T] shows time traces of Urms at Re — 79, 88, 101, 128, 
the quasi-stationary states being labeled B, C, D, E. In 
each case, the flow first reaches an axisymmetric configu- 
ration which is characterized by an initial plateau (A). Si- 
multaneously, non-axisymmetric modes with even £ and 
zonal wave number to = 2 grow fastest at an exponen- 
tial rate. Shortly before saturation occurs, the symmetry 
breaking becomes visible in physical space near the polar 
regions of the sphere as indicated in Fig. [5] (a). The previ- 
ously circular cross sections of the jets towards the poles 
are stretched elliptically as it is depicted in Fig. [5] (b-d) 
for the fiow D. The nonlinear saturation of this insta- 
bility manifests itself in a relaxation of Wims, cf. Fig. [T] 
In physical space, traveling waves with m — 2 symme- 
try emerge, their spatial structure being illustrated in 
Fig. [3l The wave propagation is directed opposite to 
the sense of rotation of the mean-flow in the individual 
hemisphere as prescribed by the driving, with the wave 
vector pointing in the zonal direction. With increasing 
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FIG. 2: Illustrations of the symmetry breaking in the flow 
towards state D. (a) Poloidal cross section of the energy den- 
sity oc during the axisymmetric phase A. Arrows indicate 
the region where elliptical stretching occurs, (b-d) Time se- 
ries depicting the spatial structure of large amplitudes of 
inside the sphere r < 0.75 as the instability grows. 
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FIG. 3: MoUweide projection of the radial component of a 
velocity fleld snapshot at 90% of the sphere's radius during the 
statistically stationary state D. Arrows indicate the direction 
of wave propagation. 



Re, the latitudinal extent of the wave grows towards the 
equator, coupling the two hemispheres and causing os- 
cillatory fluid motions (cf. Fig. [H states C and D com- 
pared to state B). As a result, the wave frequency de- 
creases with increasing Re. With Rm « 44 (see below), 
the frequencies in the lab frame turn out to be fast on 
the resistive timescale, namely /c ~ 19.08t^^ — 0.43x^7^ 
and /d ~ 16.71T-1 = 0.39r-i. For Re > 125, turbulence 
develops, and the m = 2 wave feature loses its clear struc- 
ture (state E). 

Fig. m displays time averaged energy spectra of the 
flows C and D in terms of spherical harmonic degree £ 
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FIG. 4: Time-averaged spectra of the kinetic energy density 
in the flows C and D as functions of spherical harmonic degree 
£ (a) and order m (b). 



and order to. The time averaging was performed for sev- 
eral hundred to assure convergence. The ^-spectra in 
Fig. m (a) peak at £ = 2 for both flows which is due 
to the imposed forcing. The spectra are overall similar 
except for the strength of the I = 1 modes which refer 
to global circulation, however, these modes are dynam- 
ically not important in flow C and even less in D due 
to their small amplitudes. The spectra of flow C drop 
off somewhat faster than those of flow D which is to be 
expected due to the larger value of the kinematic vis- 
cosity in run C. The spectra in terms of the order m 
shown in Fig.|4](b), however, yield differences between C 
and D concerning the amplitude of the modes with odd 
wavenumbers. In both cases, the spectra peak at to = 
which is due to the axisymmetric background flow. The 
dominant peak connected with non-axisymmetric modes 
is found at m = 2 and caused by the wave feature, fol- 
lowed by peaks at higher harmonics with even to. The 
origin of the preference for even modes in flow D remains 
unclear. During the initial bifurcation sequence, the am- 
plitudes of the odd m-modes drop off suddenly around 
Ti, « 70 after a phase spectrally similar to C before the 
statistically stationary state D (Re = 101) is reached. 
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The corresponding change of the flow becomes appar- 
ent in Fig. [T] as a change in the oscillatory behaviour 
of the respective time trace. Conversely, the flow C at 
Re = 88 starts off with odd m-modes at small amplitude 
after the initial transition A. The respective time trace 
in Fig. [1] shows regular oscillations until, at Ti, k, 210, 
odd m-modes have grown sufficiently strong to become 
dynamically important. The change in the flow towards 
the statistically stationary state C manifests itself in a 
more irregular time trace of Wrms- Flows with metastable 
states switching between the characteristics of flow C and 
D were not observed, although this possibility may exist. 
It should be noted that, when going to higher Re, the 
m-spectra tend to be smoothed. A local peak at to = 2 
is still found, and can be interpreted as the signature of 
the most unstable modes excited by the hydrodynamic 
instability. 

In the case of flow D, the temporal evolution of the 
dominant components of the spectrum is shown in the up- 
per panel of Fig. [9] The time trace of the (£ = 2, to = 0) 
axisymmetric background flow [which will be abbreviated 
(2,0) in the following] is modulated sinusoidally at fre- 
quency fo via nonlinear interactions with other modes. 
The TO = 2 wave feature corresponds to modes oscillating 
at a time-periodic non-sinusoidal pattern at an amplitude 
which is about 10% of the amplitude of the (2,0) back- 
ground flow. In the case of flow C, the picture is overall 
similar, although the time traces show different charac- 
teristics concerning their periodicity and shape. 

Additional numerical experiments were performed to 
better understand the origin of the wave phenomenon. 
Outer shear-layers which exist due to the zero-slip condi- 
tion on the velocity field prove to be crucial for the wave 
generation. This boundary condition requires the three 
vector components of the velocity field to be zero at the 
outer boundary. It is intended to be an approximation of 
the experimental reality where liquid sodium is confined 
by a solid wall. Applying the stress-free boundary condi- 
tion instead in a numerical experiment while leaving the 
forcing F unchanged fundamentally alters the nature of 
the hydrodynamic instability. The stress free boundary 
condition requires the radial component of the velocity 
field to be zero at the boundary, whereas the angular 
components evolve under the constraint of zero stress. 
As a result, a flow pattern with a strong zonal flow next 
to the boundary develops in each hemisphere. Stationary 
states with waves comparable to the flows C and D are 
not observed. These findings support the potential rel- 
evance of the results presented in this paper to existing 
dynamo experiments where boundary layers with strong 
shear exist in the flows next to the outer walls. 

Moreover, the fact that s2t2 flows were studied in the 
past by several groups, each of them using slightly differ- 
ent forcing functions or prescribed flows, motivates us to 
investigate the role of the forcing term on the wave gen- 
eration. Applying the zero-slip condition, wavy states 
similar to C and D are found in the sub-turbulent regime 
when the flow is driven by the localized body force given 



in [T^ [T3| , and, similarly, when the global s2t2 profiles 
discussed in [l^ are employed to drive the flow. In addi- 
tion to the axisymmetric component, non- axisymmetric 
TO = 2 modes are found to dominate in these flows, how- 
ever, modes with m — I and to = 3 wave numbers are 
dynamically important as well. A main difference be- 
tween the three forcing schemes consists in the distance 
between the active region where momentum is injected, 
and the outer boundary. This distance is largest in our 
impeller model, Eq. [31 and is chosen to be smaller in the 
latter forcing schemes. The decrease affects the thick- 
ness of the boundary layer and thereby the rate of strain 
which in turn may govern which zonal mode to is most 
unstable. A final comment on the influence of the forcing 
is made on the structure of the magnetic field in case of a 
dynamo. It turns out that a flow forced by Eq.[3]favors an 
TO = 1 dipole mode in the laminar case, in the intermedi- 
ate sub-turbulent range on which we focus in the present 
paper, as well as in the presence of developed turbulence. 
An TO = axial dipole as it was reported in [l^ [11] has 
not been found using Eq. [31 We have, however, numeri- 
cally conflrmed its existence in turbulent flows when the 
forcing reported in [l^ [13] is applied. 

A more detailed investigation of the hydrodynamics 
of the s2t2 flow is certainly required to shed light on 
the plausible conjecture that the previously discussed hy- 
drodynamic instability and wave formation are universal 
features of the spherically bounded s2t2 flow. The cor- 
responding study, however, goes beyond the scope of the 
present manuscript. We now proceed towards the actual 
goal of the paper, namely the investigation of the dynamo 
properties of the sub-turbulent flows C and D. 



IV. MHD INVESTIGATIONS 

The s2t2 class of dynamos has been the focus of several 
numerical inve stig ations, pioneered by the work of Dud- 
ley and James [18| who solved the induction equation for 
a stationary prescribed two-vortex flow. Later — mainly 
concentrating on the role of turbulence on the dynamo 
process — fully time dependent flows in nonlinear MHD 
simulations were studied [H, [3, [l^. Until now, how- 
ever, the effect of smoothly fluctuating fluid motions be- 
low the transition to developed turbulence has not been 
considered explicitly. This motivates us to investigate the 
temporal evolution of seed magnetic flelds in the flows C 
and D. To this end, the magnetic fleld B is initialized 
with pseudo-random noise at small amplitude, and the 
induction equation is solved numerically using four dif- 
ferent types of velocity fields v based on C and D: (i) 
time dependent flows v(x, i), i.e. the coupled Eqs. ([I])-© 
are integrated; (ii) snapshots v(x, i^), i = 1 . . .rt of the 
flows kept flxed in time; (iii) time averaged velocity fields 
V = y J^"'^'^ v{x,t)dt, with T > lO^Tj,; (iv) consecutive 
snapshots of the flows iterated at variable rates which al- 
lows to examine the effect of the wave frequency on the 
magnetic field's growth rate. Practically, the runs of type 
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FIG. 5: Energies of the magnetic field over time, from simu- 
lations (i)-(iii) based on the flow C (Re — 88, Rm — 44). 



(iv) are performed by subsequently reading previously 
saved snapshots from hard disk, while solving Eq. ([1]). 
Unless stated otherwise, all MHD runs discussed in the 
following are performed at Rm « 44. 

We first consider the flow C for which Fig. [5] summa- 
rizes the results of the numerical experiments (i)-(iii). 
Eight randomly selected snapshots — as indicated by the 
inset — show decaying magnetic fields. Moreover, the 
mean-flow is not a dynamo. Only the time dependent 
flow (i) shows a growing magnetic field solution. Clearly, 
the presence of time dependence in the velocity field in 
the form of periodic wave motion causes the magnetic 
field to grow. In fact, in the regime around Re « 100, 
the dynamo threshold of the mean flow is Rmc(v) w 55, 
which is significantly larger than the threshold of the 
time dependent fiows, being Rmc « 32. Fig. [6] displays 
the growth rate as a function of the relative wave fre- 
quency ///c obtained by performing experiment (iv). 
For comparison, the growth rate calculated from the self- 
consistent run (i) is included. A window with positive 
growth rates (///c ~ 0.55 . . .6.0) exists which contains 
two local maxima at ///c ~ 0.9 and ///c ~ 4.4 being 
symptomatic of resonant behavior. The first maximum 
happens to be in good agreement with the frequency of 
the hydrodynamic waves in run (i). Negative growth 
rates occur when the wave is either too fast or too slow. 
The latter finding refiects the fact that snapshots of the 
flow are not dynamos. Iterating through the velocity fleld 
snapshots in reversed order, i.e. inverting the direction 
of wave propagation without changing the shape of the 
waves, causes seed magnetic fields to decay in all inves- 
tigated cases in the range ///c = —10.0. . . — 0.1. The 
observation of magnetic field amplification in a spheri- 
cally bounded MHD dynamo model induced by hydro- 
dynamic waves represents the key result of our work. A 
power spectrum of the growing magnetic field obtained 
from the nonlinear MHD simulation of fiow C is given 
in Fig. [T21 The magnetic field is clearly dominated by 
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FIG. 6; Growth rate of the magnetic field, normalized by the 
negative free dipole decay rate, as a function of the relative 
wave frequency f/fc- The simulations were performed based 
on flow C at Rm = 44. 



the (1,1) dipole mode. The dynamo based on experi- 
ment (i) saturates to a stationary state without any time 
dependence in the velocity and magnetic fields. The ki- 
netic and magnetic energies are roughly in equipartition. 
Before we turn towards a discussion of the physics gov- 
erning the magnetic field growth in flow C, we perform 
an analogous investigation of the flow D. 

An intriguing phenomenon is found when the evolution 
of magnetic fields in flow D is studied. Performing run 
(i), a seed magnetic field first grows exponentially at a 
normalized rate 7d ~ 0.46, as it is indicated by the time 
trace of the magnetic energy in Fig. [7l Nonlinear feed- 
back via the Lorentz force acts starting from « 4.0 for 
a period of ~ 0.5Tcr during which the energy ratio peaks 
at E^ag/Ektn ~ 0.04 (labeled "NL" in Fig. [71). A transi- 
tion occurs to a metastable state D' which has a lifetime 
of about 2t„ w 90t^. The MHD system no longer sus- 
tains the dynamo process, as the exponential decay of 
the magnetic field at rate 7d' ~ —0.65 indicates, until it 
recovers back to the initial state D which again allows for 
magnetic field growth at 70, thus closing the cycle. In 
the following and unless stated otherwise, the labels D 
and D' refer to the growth and decay phases as indicated 
in Fig. El 

To investigate the nature of the cyclic behaviour, we 
have repeated the numerical experiments (ii)-(iv) using 
flows during both the growth phase D and the decay 
phase D'. The mean flows and field snapshots were ob- 
tained from run (i) at times sufficiently far away from the 
transitions and while the magnetic field was negligibly 
small. Experiment (ii) shows growing and decaying solu- 
tions for the magnetic field in both cases, D and D'. This 
result does not disagree with the result obtained from 
flow C at Rm = 44 since a similar situation can easily 
be constructed using D by decreasing Rm below 44, i.e. 
all selected snapshots in experiment (ii) are not dynamos 
whereas a random seed magnetic field grows in time. The 
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FIG. 7: Time traces of kinetic and magnetic energy, the flow 
initiaUy being in state D (Rm = 44). The inset shows the 
kinetic energy during the transition from growth to decay. 



mean-flows during phase D and D' are identical by visual 
inspection, which explains the finding that experiment 
(iii) yields identical decay rates 7dm ~ —0.31. Note that 
during the phase D' in run (i), the magnetic field de- 
cays twice as fast as it does in the mean flow. Finally, 
experiment (iv) yields resonance behaviour qualitatively 
similar to flow C. To shed more light on the nature of 
the transitions, further investigations are presented in 
the following. The main question to be adressed is how 
the velocity field, the magnetic field, or the combination 
of both is altered during the transitions. 

"Self-killing" (though — to our knowledge — no "self- 
recovering") nonlinear dynamos have been reported pre- 
viously [20, [2l| ■ Their mechanism is based on a modifi- 
cation of the flow due to the action of the Lorentz force 
away from an initial state which has supported magnetic 
field growth towards a second solution of the nonlin- 
ear Navier-Stokes-Eq. Q which prevents magnetic field 
growth. During the transition from D to D' additional 
exponentially decaying smooth fluid motions on a slow 
timescale are excited, as the inset in Fig. [7] indicates. 
However, the magnetic field continues to decay at con- 
stant rate 7d' while the fluid perturbation is damped to 
very smafl amplitude, indicating that this perturbation 
might be unimportant. This hypothesis will be confirmed 
by the investigations presented below. 

To scrutinize the transition, we have performed the fol- 
lowing numerical experiments (a)-(c) using the solutions 
vnl and/or Bnl of the system at time t„ « 4.2 as initial 
conditions in the MHD equations. The field Bnl was 
rescaled to infinitesimal amplitude, comparable to the 
seed fields used in previous runs. Hence, the structure 
of the magnetic field is preserved whereas an immediate 
backreaction on the flow is avoided. The experiments 
and results are as follows, (a) Integration of the MHD 
equations using vnl and the rescaled Bnl as initial con- 
ditions exhibits the decay phase at rate 7d' [which is of 



the same duration as in (i)], followed by the transition 
to growth, (b) Integration starting from vnl and an in- 
finitesimal pseudo-random magnetic field yields magnetic 
field growth at rate 70. (c) Integration starting from a 
velocity field obtained during D' different from vnl and 
from the rescaled Bnl (i-e. introducing a phase shift be- 
tween the self-consistently adjusted configurations of the 
fields) leads to the same result as in run (b). 

We have further investigated the transitions by exam- 
ining explicitly the temporal characteristics of the spec- 
trum. To this end. Fig. [8] displays time traces of the sum 
of the volume-integrated poloidal and toroidal energies 
in the dominant modes during the transition from D to 
D'. As it was already pointed out in the previous sec- 
tion on hydrodynamics, the (2,0) time trace of flow D 
is modulated sinusoidally. Hence, it is suited to derive 
the phase relative to a mode which oscillates similarly 
in time. In Fig. [5] (a), the dominant modes of kinetic 
energy obtained from the MHD simulation are plotted 
over time. These are the (2,0) mode which is responsi- 
ble for the axisymmetric background flow, and the three 
strongest (m = 2) modes which account for the bulk 
of the non-axisymmetric wave feature. In addition, the 
(2,0) time trace from a simulation without magnetic field 
is shown. The generating run was started using an initial 
condition obtained from the MHD simulation at t„ = 2.0 
when the magnetic field was negligibly small. It is seen 
easily by comparing the (2,0) curves that the nonlinear 
feedback due to the Lorentz force causes a phase shift 
by TT relative to the unperturbed flow. Except for this 
phase shift, the spectral characteristics of the flow re- 
main unaffected. Fig. [5] (b) shows the dominant mag- 
netic energy modes during the same time interval, and in 
addition, the time trace of the (2,0) kinetic energy mode 
from Fig. [5] (a). Before the transition takes place, the 
oscillations in the dominant modes of the magnetic and 
the velocity fleld are in phase, as it is indicated by ver- 
tical lines. After the transition, the curves oscillate out 
of phase by tt. In particular, the magnetic energy modes 
continue to oscillate at their characteristic frequencies 
during the transition whereas the kinetic energy modes 
experience a phase shift by tt. This phase shift is related 
to a change in the relative alignment of the velocity and 
magnetic fields which will be discussed in detail in the 
following paragraph. The relevant spectral components 
during the transition from D' to D are displayed in Fig.[51 
It is readily seen in subfigure (a) that the velocity field 
modes are completely unaffected by the change from de- 
cay to growth of the magnetic field. The magnetic time 
traces switch from out-of-phase oscillations to oscillations 
which are in phase with the dominant hydrodynamic os- 
cillations as subfigure (b) indicates. 

It is illustrative to translate the phase shifts to the 
respective changes in physical space. To this end, we 
begin with an investigation of the velocity field during 
the transition NL. The panels (a) and (c) in Fig. [TUl show 
snapshots of the radial component of the velocity field 
taken at maxima of the (2,0) kinetic energy time trace. 
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FIG. 8: (a) Volume-integrated spectral energies as functions 
of time of the 4 dominant modes in flow D during the tran- 
sition from D to D' in the full MHD simulation, c.f. Fig. [7] 
For comparison, the time trace of the (2,0) mode of a mag- 
netically unperturbed, purely hydrodynamic simulation of the 
same flow is included, rescaled to 50% of its amplitude. The 
scaling is indicated by an asterisk (*). Vertical bars indicate 
the phase shift which is imposed on the hydrodynamic os- 
cillation, (b) Spectral energies as functions of time of the 3 
dominant modes of the magnetic fleld during the transition 
from D to D'. For comparison, the (2,0) mode of the veloc- 
ity fleld is displayed, cf. subfigure (a). Vertical bars illustrate 
a phase shift which develops between the oscillations of the 
magnetic energy and the oscillations of the kinetic energy. 
The magnetic energy (*) was scaled by a factor of 10. 

before and after the transition, viz. Fig. [8] (a). For the fol- 
lowing discussion, it is useful to introduce a virtual plane 
on a circle of longitude which is defined by the z-axis 
and the longitudinal coordinate at which the counter- 
propagating wave features in the upper and lower hemi- 
spheres "meet." The position of the plane varies between 
runs due to the rotational symmetry the flow has about 
the z-axis before it becomes hydrodynamically unstable. 
However, once states like C or D are reached, this plane 
stays fixed in space since the waves counter-propagate 
at the same velocity in each hemisphere. Relating this 
virtual plane to the oscillations in Fig. [5] (a), the wave fea- 



FIG. 9: (a) Dominant modes of the velocity field during the 
transition from D' to D, viz. Fig. [H](a). (b) Dominant modes 
of the magnetic field during the transition from D' to D, viz. 
Fig. |5] (b) . The (2,0) velocity mode from subfigure (a) is 
displayed for the purpose of illustrating the phase angle rela- 
tive to the magnetic dipole mode. For practical reasons, the 
amplitude of the kinetic (magnetic) energy was scaled by a 
factor of 9"^ (10^°). 



tures reach it at maxima of the (2,0) kinetic energy time 
trace. From Figs. [10] (a) and (c) it is evident that the 
nonlinear feedback during NL due to the Lorentz force 
rotates the plane by 7r/2 about the z-axis which trans- 
lates to the phase shift by tt relative to the unperturbed 
flow as indicated in Fig. [5] (a). The panels (b) and (d) of 
Fig. [To] show the radial magnetic field snapshots associ- 
ated with the velocity fields shown in the panels (a) and 
(c). The magnetic field is dominated by an (1,1) equato- 
rial dipole mode which is evident from the fiux patches 
peaking on the equator. The magnetic field's geometry 
changes periodically in time during wave motion, how- 
ever, the equatorial dipole stays fixed sufficiently far from 
the transitions which can be understood from the same 
symmetry argument as in the case of the aforementioned 
plane. The magnetic equatorial dipole does not change 
its position during NL as a comparison of the panels (b) 
and (d) in Fig. [TO] yields. Comparing the panels (a) and 
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FIG. 10: Radial vector component of the velocity field (a, c) 
and the normalized magnetic field (b, d) at 90% of the sphere's 
radius before (a, b) and after (c, d) the first transition from 
growth to decay, labeled NL. The snapshots were taken at 
maxima of the (2,0) kinetic energy time trace, viz. Fig.[8](b). 
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FIG. 11: Radial vector component of the normalized magnetic 
field at 90% of the sphere's radius during the second growth 
phase (a), the second decay phase (b), and the third growth 
phase (c). 



(b) in Fig. 1101 the plane characterizing the wave motion 
is co-aligned with the equatorial dipole before the tran- 
sition. After NL, the plane and the dipole are oriented 
perpendicular to each other due to the changes in the 
velocity field, viz. Fig. [TO] (c) and (d). Obviously, the 
ahgnment of the velocity field relative to the magnetic 
field is decisive if the magnetic field grows or decays in 
time. Let us continue with following the changes to the 
fields during subsequent transitions. Fig. II II shows snap- 
shots of the normalized radial magnetic fields during the 
second growth phase, the second decay phase, and the 
third growth phase. During the transition from decay to 
growth at « 6.7, the equatorial magnetic dipole ro- 
tates by 7r/2 as it is evident from comparing Figs. [TO] (d) 
and nil (a), whereas the flow remains in the state indi- 



cated by Figs. [TO] (c) and [9] (a). Hence, the magnetic 
dipole component and the plane characterizing the wave 
motion are parallel to each other again which results in 
magnetic field growth. Once the magnetic field has grown 
sufficiently strong which is the case at « 10.1 the wave 
motion is pushed towards the initial state as depicted in 
Fig. [TO] (a) while the equatorial dipole keeps its position, 
cf. Fig.[TT](b). The magnetic field decays and starts grow- 
ing again when its equatorial dipole component rotates 
by tt/2 to reach co-alignment with the plane character- 
izing the wave motion, cf. Figs. [TO] (a) and [TT] (c). In 
summary, the equatorial dipole tilts by 7r/2 about the 
z-axis between subsequent growth phases which adds up 
to a field reversal after two cycles, cf. Figs. [TO] (b) and 
[TT](c). The magnetic induction Eq. [fallows for such field 
reversals since if B is a solution for a given flow, — B is 
a valid solution as well. Whenever the magnetic field re- 
acts back on the velocity field, the flow is forced to switch 
between the configurations shown in Figs. [TO] (a) and (c). 
The dynamo based on flow C operates similarly to flow 
D during the growth phase as it is shown in Figs. [TO] (a) 
and (b), however, does not show cyclic behaviour. 

The power spectrum of the growing magnetic fleld in 
system D is given in Fig. 1121 The power spectrum of 
the decaying magnetic field is virtually undistinguishable 
and was therefore not included. Compared to C, the 
odd m modes are more pronounced which is due to the 
preference for even m modes in flow D. 

Nonlinear saturation to a stationary state with 
Emag/ Ekin ~ 1-2 is reached above Rm = 51. Time 
traces of the dominant kinetic energy modes during satu- 
ration are shown in Fig. 1131 Compared to the previously 
discussed cyclic dynamo the system saturates towards a 
completely stationary state without any oscillations sim- 
ilar to the saturated dynamo state of flow C at Rm = 44. 
We now turn to a discussion of the physics results. 



V. DISCUSSION 

Our numerical investigations indicate that magnetic 
field amplification generated by time periodic non- 
axisymmetric m = 2 hydrodynamic waves is relevant to 
the large-scale dynamo process inside a spherical fluid 
conductor. These waves emerge due to a hydrodynamic 
instability which breaks the rotational symmetry of the 
flow. The flow itself is driven by an axisymmetric body 
force and subject to the zero-slip outer boundary condi- 
tion in order to mimic a situation similar to the Madison 
Dynamo Experiment. Several bifurcation sequences to- 
wards different statistically stationary flows with m — 2 
waves are observed depending on the fluid Reynolds num- 
ber. The numerical investigations are performed in the 
intermediate regime between the laminar flow and devel- 
oped turbulence which facilitates the identiflcation of the 
dominant physical effects. 

Let us briefly reiterate the key results which were ob- 
tained by solving the induction equation using the afore- 
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Order m 

FIG. 12: Normalized and time-averaged spectra of the mag- 
netic energy density during the kinematic (growth) phase in 
the flows C and D as functions of spherical harmonic degree 
i (a) and order m (b). 
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FIG. 13: Saturation of a dynamo based on flow D at Rm = 51. 
The system evolves towards a completely stationary state as 
it is illustrated by time traces of the dominant kinetic (a) and 
magnetic (b) modes. 



mentioned sub-turbulent flows. At a given magnetic 
Reynolds number, the time-averaged flows as well as 
snapshots of the flows are not dynamos whereas the time- 
dependent flows with wave motion do show dynamo ac- 
tion. By variation of the wave frequency, maxima in the 
efficiency of the dynamo process are identified, as well as 
the breakdown of the dynamo for frequencies which are 
too low or too high. Moreover, one of the several possi- 
ble statistically stationary states of the flow shows cyclic 
dynamo action. We find that the transitions between 
magnetic field growth and decay crucially depend on the 
phase angle between characteristic oscillations in the ve- 
locity and magnetic energy modes which translates to 
changes in the relative alignment of the fields in physical 
space. Hence, system D is different from the self- killing 
dynamos which have been found previously [13, |2l[ . 

These phenomena are symptomatic of dynamo action 
based on non- normal growth p^ . Non-normal growth 
refers to perpetual growth due to the mixing of non- 
orthogonal eigenstates of the time-dependent linear op- 
erator Cii) of the induction equation, which can formally 
be written as d'Qjdi — >C(t)B. If the time dependence in 



C(t) given by a non-stationary fiow is adequate, perpet- 
ual growth of magnetic field modes is possible even if the 
magnetic field would decay in the absence of time depen- 
dence. Such mixing of non-orthogonal eigenstates was 
shown to act as an efficient driver for dynamo action in 
Ref. 16], thereby identifying a novel dynamo mechanism 
different from Lagrangian chaos or the a-effect. Using an 
infinite 2D periodic drifting velocity field which is no dy- 
namo without drift at a given Rm, it was demonstrated 
that the velocity field supports magnetic field growth 
for drift velocities not being too fast or too slow. The 
analogy between our results and the results presented 
in Ref. [l^ strongly supports the conjecture that non- 
normal growth is responsible for the magnetic field gen- 
eration in the sifl flow with wave motion. Interestingly, 
flow C exhibits a wave frequency close to a maximum in 
dynamo efficiency, cf. Fig. [6l This happens by chance 
since the wave feature is a hydrodynamic phenomenon 
and, hence, is independent of a weak seed magnetic field. 

The observed saturation of the s2i2 dynamo to a state 
without wave motion confirms the relevance of the ad- 
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justment of the wave frequency during the saturation 
process ^] , cf. Fig. [121 

Increasing the Reynolds number above 125 leads to a 
superposition and interaction of the m — 2 waves with 
structures at larger wavenumbers. Turbulence develops 
which manifests itself in an increase of the critical mag- 
netic Reynolds number Rmc (Re), as it was shown numer- 
ically for s2t2 systems l^, 13, [Si ■ This increase in the 
stability curve due to turbulence is a big challenge for 
laboratory dynamo experiments which use simply con- 
nected flows H, 0, [E S 0, 111 ■ Liquid sodium which is pre- 
ferred as a working fluid has a magnetic Prandtl number 
Pm = Rm/Re of the order of 10~^ [T^. Hence, in order 
to achieve magnetic Reynolds numbers sufficiently large 
for self-excitation, the experimental flows are vigorously 
turbulent with Reynolds numbers Re > 10^, inaccessible 
to any numerical simulation. Therefore, only extrapo- 
lations of Rmc(Re) can be done. In numerical studies 
of dynamos with Pm < 1 conducted in infinite periodic 
boxes, which are usually computationally less expensive 
than simulations in spherical geometry but lack realistic 
boundary conditions, the stability curve Rmc(Re) was 
shown to exhibit a saturation to a plateau after an ini- 
tial increase for large-scale dynamos d, [13, [O, Ull , and 
similarly for the fundamentally different small-scale dy- 
namo driven by non-helical random forcing, cf. [ 2^ and 
the references therein. In spherically bounded large-scale 
dynamos like the system we are focusing on in the present 
study, the question if and at which level the dynamo 
threshold saturates to a plateau at large Re has yet to 
be investigated. Our study shows that field generation 
due to wave motion is important in the context of the 
dynamo instability. Beneficial effects are shown to ex- 
ist in the sub-turbulent regime at low Reynolds numbers 
of about 100. The question most interesting to experi- 
menters is, of course, if these results can be transferred 
to laboratory dynamos. Simulation results obtained for 
Reynolds numbers above 1000 are encouraging since we 
still find a pronounced peak at m = 2 in the spectra of 
turbulent velocity fields. This indicates that remnants 
of the m = 2 waves could be present in the form of co- 
herent structures in experimental flows and potentially 
act as a driver for magnetic field growth. A modification 
one might suggest in order to exploit the findings of this 
paper to improve s2t2 setups is the installation of de- 
fiector baffles in front of the impellers, designed to force 
flow components with m = 2 symmetry on top of the 
axisymmetric background flow. These baffles could then 
be rotated about the impeller axis at variable frequency 
which would allow to search for an optimum in frequency 
similar to the study presented in Fig. [6l 

Growing and transiently decaying magnetic solutions 
are supported by the flow D. The question if the mag- 
netic field grows is related to an alignment problem of 
the velocity and magnetic fields. The configuration of v 
and B depicted in Figs.[Tn](c) and (d) constitutes a poor 



state vector of the system since subsequently applying 
the linear operator C{t) of the induction equation on B 
in order to evolve the system forward in time — which is 
done by the numerical code — causes the magnetic field 
to decay. Note that the time-dependence of C{t) is given 
by the temporal evolution of the velocity field starting 
from the configuration shown in Fig. [TO] (c) . The modi- 
fication of the flow during the transition from growth to 
decay is certainly caused by the Lorentz force, however, 
the reason for the rotation of the magnetic dipole during 
the transition from decay to growth is not obvious. Slow 
diffusion of the magnetic field into the preferred position 
which is then amplified to finally take over is a possible 
explanation. 

Let us finally discuss the relation of the cyclic dynamo 
D to the competing axial and transverse dipole modes 
which are reported for an s2t2 flow in 14]. Ref. [l3| de- 
scribes the generation of an axisymmetric magnetic field 
component via non-axisymmetric turbulent fluctuations 
in the presence of an axisymmetric s2t2 background flow. 
It is shown that, first, an equatorial dipole magnetic 
field is generated which creates via Lorentz braking an 
m = 2 velocity field component which in turn generates 
the TO = magnetic field. In the following, the to = 
and TO = 1 magnetic modes compete via a modulation 
of the amplitude of the m — 2 velocity mode. However, 
these competing modes are fundamentally different from 
the cyclic dynamo D which is presented in the paper at 
hand. First, the to = 2 non-axisymmetric hydrodynamic 
waves in our model emerge independently of a magnetic 
field. Second, when a weak seed magnetic field is in- 
troduced, the flows with to = 2 wave motion generate 
a dominant to = 1 transverse dipole component. The 
axisymmetric components of the magnetic field are rel- 
atively weak as the power spectrum in Fig. [12] (b) in- 
dicates. Third, there are no competing modes observed. 
The relative alignment of the magnetic and velocity fields 
is decisive if the magnetic field grows or decays. In sum- 
mary, our results are by no means contradictory to the 
results reported in [l3| • Differences in the forcing scheme 
and in the turbulence intensity exist which are crucial if 
an axisymmetric dipole is generated. This is pointed out 
in section Hill 

In conclusion, magnetic field generation due to wave 
motion was shown to exist in a spherical s2t2 dynamo 
model. The transfer of these results to experimental 
fiows is not straightforward due to turbulence, neverthe- 
less wave effects could play a role in laboratory dynamos 
and should be kept in mind when modifying existing or 
designing new experiments. 
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